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5-H , Abstract. In this paper, we study the discrete Morse flow for the Ricci flow on 

football, which is the 2-sphere with removed north and south poles and with the 
metric ^o of constant scalar curvature, and and for Porous media equation on a 
bounded regular domain in the plane. We show that with a suitable assumption 
about ^(0) we have a weak approximated discrete Morse flow for the approxi- 
mated Ricci flow and Porous media equation on any time intervals. 
Mathematics Subject Classification 2000: 53Cxx,35Jxx 
r^ ' Keywords: discrete Morse flow, Ricci flow. Porous-media equation, conical 

>-s , singularities 

^' 
-)— » 

1. Introduction 

There are relative few results about computational models for the Ricci flow 
in two dimensions. The purpose of this paper is to try this area by giving some 
^ I approximated computational models, namely the discrete Morse flow for the 2- 

^ ■ d Ricci flow. We shall first consider the Porous-media equation on a bounded 

CN ! regular domain in the plane. As is well-known that the limiting equation of the 

^ ' Porous-Media equation is the Ricci flow on the domain. Then we may consider 

cn . this discrete flow as the approximated computational scheme for the Ricci flow. 

^ I We also consider two diff"erent modelings of the Ricci flow on the singular surface, 

the American football. Our method can also be worked out in a similar way in the 
regular spherical surfaces and other surfaces. 

Since the geometric and analytic parts of singular surfaces are not well-known, 
rS ■ we now recall the geometry of a special singular surface, namely, the American 

j^ I football. Let S = Sq he the sphere with north and south poles removed and with 



the metric 



gQ = dr^ + (a sin rf-dG^, 



where a € (0, 1), < r < ;r such that r = corresponding the north pole, and 
< 9 < 2n. With this metric, we can see that ^o has two singularities of equal 
angle lira at the north and south poles. By a direct computation we know that 
the scalar curvature of the metric go is equal to 2. Recall here that the Laplacian 
operator of go on function is 

A a2 cosr 1 2 

Ao - 5; + - — dr + - — -—^di 
sm r {a sm ry 
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and the area element is 

dA = a sin rdrdO. 
The area of iS,go) is \A\ = Ana. The aim of this note is to study the normalized 
Ricci flow on (5,^0): 

dtg = ip- R)g 
where g = g{t) = e^"^'^go, R = R{g) is the scalar curvature of the metric g, and p 
is some real constant (and we may choose it such that the area of the flow g{t) is 
constant). We may write the evolution equation of the Ricci flow as 

e"dte" = ^pe^" + Am - 1, 

where A is the Laplacian operator of the metric ^o- The standard example is that 
for the football metric go = dr^ + (asmr)^d9^ (0 < o- < 1) the (un-normalized) 
Ricci flow is g{t) = (1 - 2t)go for ? < ^, which extinct at ? = ^. 

For the construction of the discrete Morse flow for the 2-d Ricci flow, we shall 
rely on some variational structure related. Hence we need some analytical part of 
the related variational functional. Here we recall one important consequence of a 
Moser type inequality proven by W.Chen and C.Li [6|. For any € H\S), 



(1) 

and the functional 



f e'>'dA < Cjexpi— — ( \V^\^dA + — \ MA] 
Js ' 167rQ'J5 \S\J 



(peH\S)^ j e'f'dA 



is continuous with respect to the weak convergence in H^ . For our purpose, we 
reformulate the inequality above as For any u e H^{S), 

(2) -f e^''dA < C2exp{ f |VmP JA + -f 2udA} 

Js Js Js 

with j- udA = jji / udA. We shall use this fact to get a weak solution to the Ricci 

flow on (5, go)- 

The precise results will be stated and proved in the following three sections. 
Similar results for Yamabe flow and higher dimensional Porous-media equations 
are also true. The plan of this paper is below. In section |2l we study the existence 
of a weak solution to the porous media equation in a bounded regular domain. 
In section [3l we introduce the method of the discrete Morse flow method to an 
approximated ricci flow in the singular surface with symmetry the initial data. We 
study a perturbated Ricci flow model in the last section. 

2. DISCRETE Morse flow for the Porous-media equation 

The existence of the weak solution to the Porous-media equation can be done 
by using the Galerkin method [2]. Other related methods can be found in [8]. Here 
we propose a new method, which is the adapted discrete Morse flow method lfT3l . 

To make our idea more clear, we start from the discrete Morse flow for the 
Porous-Media equation on the bounded regular domain Q. in the plane R^. In some 
sense our domain be be any one with the Sobolev imbedding theorem holds true. 
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Assume that m > 0. Given an initial regular data mq and any T > 0. We consider 
the porous media equation 

(3) ^tu = ^u!", Qx[0,r], 

with the initial data u\t=Q = uq and with the boundary condition u{t) = uq on dO. x t. 
This equation has a very close relation with the Ricci flow. In fact, by taking the 
limit m — > 0, the limiting equation of (O is 

(4) dfU - A log u, 

which is the Ricci flow on the plane R^. 

To introduce the discrete Morse flow for (|3) we set v = u'" and let a = l/m. 
Then ([3]) is induced into 

(5) drv" = Av 

with the initial data v = vq = u" at t = 0. We assume that vq € H^{Q.). The key 
idea is that we set a = 2/3 - 1 with /? > 1 /2 and make 

where C^ > such that 2^0? = 1- 

In below, we assume that /? > 1 and make the notation that v^ = Ivl'' for any v. 
For any N > Ibea. large integer and for any T > 0,let 

h ^ TIN, t„ =nh, n = 0, 1,2, ...,A^. 

Assume that we have constructed Vj e H^(S), < j < n - I and m„_i is a 
minimizer of the functional 

In-iiu) = ^ J 1/ - ^^^_/dx + i J \Vv\^dx 
on H--{ve H\Q.y,v - vq e H^iQ.)]. Define 

Iniu) - ^ J 1/ - /„^|2jx +^J Wvfdx 

on H. It is clear that the infimum is finite and by using the Poincare inequality to v- 
vo, any minimizing sequence is bounded in H. We remark here that we may choose 
the functions in the minimizing sequence non-negative since our functional is even 
functional HI. By the direct method, we know that /„ has an unique minimizer m„ 
in H, which satisfies that 

h ^ "-'^ 

with the uniform energy bound 



.i\^dx<C. 



We inductively define v^it) e H^ for t e [-h, T] such that for n = I, ...,N, 

VNit) = V„, 
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on [tn-i , t,,] and VN{t) - vq on [-h, 0]. We define, for \ <n < N, 

and VA?(f) = VN{t-h). Taking tlie convergent subsequence in the weak star topology 
in H, we know that the limit v £ // satisfies 

Cp^P-^^,^P = Av 

in the distributional sense. The latter is the weak form of the equation ([5]). 
Then we have proven the following result. 

Theorem 1. Assume that Q. is a regular domain in the plane R^. Given any T > 
and m > 0. Assume that the initial data vq = uJ'" € H^{D.). Let a = l/m = 2/3 — I 
with p > 1. Then there is at least one weak solution v e L°°H^{£l) to the porous 
media equation dJ]). 

Recall here that a mapping v : [0, T] — > H^ (Q) is said a weak solution to © if 
V ee U°H^{Q.) satisfies the evolution equation (|5]l in the distributional sense and 
lim^^o u{t) = uq in L^{Q) and u{t, x) - uq{x) for x e dQ. in the trace sense. The 
latter will be simply said that u{t) has the initial data mq- 

3. H^ WEAK SOLUTION TO RiCCI FLOW ON {S , g()) WITH SYMMETRY 

We now choose p - . , and fix any positive constant T > 0. Fix any 

Mo £ H^ with u{x) = u{-x), which is symmetric about r = n/2. Then the Ricci 
flow equation under consideration is 

(6) e"d,e" - Am - 1 + , in S x (0, T] 

j e^"dA 

with the initial data m(0) = uq. The symmetry condition is used to get the com- 
pactness required for the minimization process below. We now introduce a new 
concept of H^ weak solution to the Ricci flow Q. 

Definition 2. We say u{t) e H\ t e [0, T] is a weak solution to the Ricci flow Q 
with the initial data uq if it satisfies ([61) in the sense of distribution and with the 
bounds 

sup -f \Vuf- + ludA - log +■ e^"dA < \ |VmoP + lu^dA - log \ e^'^'dA. 

and 

f \dte"^'^\^ < ^ r IVmoP + 2uodA - ^ log f e^"°dA. 

J5x[o,r] '^ Js ^ Js 

In below, we may assume that mq is smooth (i.e., uq £ C^) and normalize it 
such that the average "o = _;£"- iji C uodvg^ , otherwise, we choose an smoothly 
approximation of uq and passing to subsequence limit. 

For any N > Ihe a. large integer and for any T > 0,let 

h = T/N, tn = nh, n = 0,1,2, ...,N. 
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1/ 



Assume that we have constructed Uj € H (S), < j < n - I and m„_i is a 
minimizer of the functional 

Jn-i{u) = ^j \e" - e""-'\^dA +^-j \^uUa - ^ log J e^'^dA 



over X - [u e H ,u -0} and with the uniform bound 



■f \Vun-\\^dA - log f e^"dA < C. 
Js Js 



We want to get another m„ with same bound by using the discrete Morse method 
due to Rothe . Define the functional 

Jniu) = ^j k" - e""-' Ua + ^j \VuUa - ^ log J e^'^dA 

in A. If Un is a minimizer of the functional on A, then it is clear that its Euler- 
Lagrange equation of this functional is 

1 „2Un 

(7) e"" -(^"« - e"'-' ) = ^Un- K + 7^-^, on S, 

n j- e^"" 

where 

Js h 
and the minimizer satisfies the following estimate 

(8) Uun) < Uun-i) = ^ r iVM„-i P JA " ^ log j^ e^""-'dA. 

We want to minimize this functional on A. First we need to know that the 
infimum is finite. This follows from the fact that after using the inequality ([T]), the 
leading term in the functional /„ is ^ C k" - e"" ' ^dA. 

To get the minimizer, we only need to show that the minimizing sequence is 
bounded in A. Right from the the relation 

inf 7„(m) < Jn{Un-l) < C, 

hHs) 
we know that the minimizing sequence (u = Unu) satisfies that 



J_ / |e" _ e««-i \^dA + \f \^u\^dA - ^ log -f 



e^"dA < C. 



Using Chen-Li's Moser type inequality we know that the minimizing sequence is 
uniformly bounded in H^. By this we can pass to subsequence to get a weakly 
convergent subsequence with its limit m„ as the minimizer. 

Once this is done, we inductively define unQ) e H^ for t e [-h, T] such that for 
n = l,...,N, 

Unit) = Un, ^Nit) - '^n, 

on [tn~i, t„] and u^it) = uq on [-h, 0]. We define, ior I < n < N, 

d,e"^(^-> = he"" - e""-'), f€[f„_i,f„], 
n 
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and mtvCO = ui\j{t - h). The estimate ([H) gives us that 



sup f \Vun\^dA - - log -f e^""dA < \ |Vmo|Va - - log f e^^oj^ 

and 

r \dte"''^Y <lf IVmoP^A - ^ log f e^"'^dA + C,,. 

JSx[0,T] '^ Js ^ Js 

The uniform upper bound of L IVma^PJA follows from the uniform bound of 
f g"iv and the assumption that u{x) = u{-x) where the reflection is about r - njl. 
With the latter symmetry assumption we have the better inequality in the sense of 
Moser that 

-f e^"dA < C2exp{- -f IVwpJA + -f ludA) 

which is proved by W.Chen ||5l (see theorem II there). It is here that the symmetry 
condition plays the role. 

Note that [un] is bounded in L~//' and \d,e"''\ is bounded in L^{\0, t]xS) for 
any ? > 0. Then we may pass to weakly convergent subsequence of u^it) in (|7]) to 
get the weak limit u{t) e L^H\S) n L\[0, t]xS) such that 

2m 

(9) e"dte" = Au- A{t) + -— , 

Is ^ "^^ 
where 

A(t) - lim AN{t) - 1 - f e'^dte". 

The latter is called the relaxation parameter of the ricci flow and the open problem 
is to prove that it decays to zero as time going to infinity. 
Hence we have proven the following assertion. 

Theorem 3. For any positive constant T > and any uq € H^ with u(x) = u(—x) 
with respect to r - n/l, there exists at least one weak solution to ^ on S x [0, T] 
in the sense of the weak form dP]). 

Our argument above also work in the case of the standard sphere. It may be 
used to define the weak solution to the Kaehler-Ricci flow. It is an open question 
if we have stability result or uniqueness result for the weak solution above. We 
remark that it is possible to get a little more regularity of the weak solution for 
smooth initial data uq. For higher dimensions, except the concept of Ricci flow 
with surgery, it is not clear if one can define the weak Ricci flow. 

4. Another discrete Morse flow 

Fix A G (0, 1). This parameter will play the regularization role in the discrete 
Ricci flow below. We study the following modified Ricci flow 

2u 

(10) e"dte" = Am - /I + , in S x (0, T] 

fe^"dA 
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with the initial data m(0) - uq € H^. In this case, the variation structure related to 
the right side of (ITOl ) is 

I(u) = - f (|Vm|2 + 2Au)dA - -log f e^"dA. 

The advantage of this function is that the inequality ([2]) implies that 



■f, 



(11) liu) > 2){A - I) -r udA 
and 

(12) I{u) > — — -f- |VmP JA — log -f- e^^dA. 

These two facts will help us to make the infimum of the functional J (see below) 
in H^ be finite. 

As in last section, for any N > Ihta large integer and for any T > 0,let 

h = TIN, tn = nh, n = 0, 1,2, ...,A^. 

Assume that we have constructed iij € H\S), < j < n - \ and m„_i is a 
minimizer of the functional 

Jn-\{u) = — f k" - e""-2pc/A + - f (|VmP + 2Au)dA - - log f e^"dA 
2h Js 2 Js 2 Js 

over H^ and with the uniform bound 

■f (|Vm„_iP + 2Au)dA - log f e^"dA < C. 
Js Js 



Let 



4|i'" 



/(„) - — f |e" - e""-' \^dA + I{u) 



on //^ Because of the inequality Q we can easily show as before that the mini- 
mizing sequence of J is bounded in H^. In fact, if the infimum is oo, then we may 
have sequence wj € H^ such that Jiwj) -^ -oo. By (fT2l) or (ITTI ) we know that 
r e^'^V _> oo or _^ Wj —^ oo, which is impossible since in this case the leading term 

is 

1 r 

„«n-l |2 



^i' 



Is 

Hence the infimum is finite. By this we can pass to subsequence to get a weakly 
convergent subsequence with its limit m„ as the minimizer of J. The Euler-lagrange 
equation of J is 

1 g2a„ 

(13) e"»-(e"« - e"'-') = Am,, -A + , on S. 

h + e^"" 

By the standard elliptic regularity theory we know that u,, is smooth. 

Then we inductively define mw(0 € H^ for t € [-h, T] such that for « = I, ...,N, 

U!\/{t) = Ufi, 
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on [tn-i , t,,] and ma'CO = "o on [-h, 0]. We define, for I < n < N, 

h 
the minimizing property of m„ gives us the desired energy bound to obtain a con- 
vergent subsequence un and the limit m e L'^H^{S) n L^([0, t]x S) such that 

(14) e"5,e" - Am - /I + 



j's ^^"^^ 



Hence we have the following result. 



Theorem 4. For any positive constant T > and any mq £ //\ ?/jere exists at least 
one weak solution to U0\) on S x [0, T] in the sense of the weak form U4\l . 



Here the weak solution of (flOl ) is defined similar to the definition |2l 
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